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The data for type la supernovae obtained by the High-z SN Search Team and Supernova 
Cosmology Project are analyzed using inhomogeneous cosmological models with a local void 
^ , on scales of about 200 Mpc, to derive the best-fit values of model parameters and the 

confidence contours. The x 2 fitting is found to be slightly better than that in homogeneous 
models. It is shown that (1) the best-fit values are most sensitive to the ratio R of the 
outer Hubble constant, Hq 1 , to the inner Hubble constant, Hq, (2) the best-fit outer density 
parameter, and cosmological constant parameter, Aq 1 , are, respectively, increasing and 
decreasing functions of R, and (3) (fij 1 ) -^o 1 ) can t> e (lj 0) f° r R ~ 0-8- Moreover, it is shown 
that these models are naturally consistent with the new supernova data (SN1997ff) with 
> ' z= 1.7. 
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§1. Introduction 



In present day cosmology, the most important observations are those of the 
[magnitude m - redshift z] relation for type la supernovae (SNIa) and the CMB 
anisotropy. In the former observations, SNIa play the role of the best standard 
Q-f candles, and two groups, the High-z SN Search TeamtP'cP'lf and the Supernova 

Cosmology Project, B have to this time observed 50 and 60 SNIa, respectively. Their 
results suggest that a significant amount of dark energy fills our universe, under the 
assumption of its homogeneity and isotropy. 

Inhomogeneous models with a local void on scales of about 200 Mpc have been 
studied by the present author EP in connection with puzzling behavior of cosmic bulk 
flows cP'tP and it was subsequently shown that it may be possible to explain the accel- 
erating behavior of high-z SNIa in [m, z] relation without a cosmological constant.^ 
A historical survey of works concerned with local voids is given in a previous work.H* 
Because these inhomogeneous models include many parameters (such as the inner 
and outer values of the Hubble constant and the density parameter, the cosmologi- 
cal constant, and the boundary radius), we recently examined the dependence of the 
above relation on these parameters, in comparison with the relations in homogeneous 
models with (fl , A ) = (0.3,0.7) and (0.3, 0.0). EP 

Observationally, recent galactic redshift surveys liiP'Ej)'t3'lliP show that in the 
region approximately 200 - 300/i _1 Mpc from us, the distribution of galaxies may be 
inhomogeneous. Moreover, a large-scale inhomogeneity suggesting the existence of a 
wall around the void on scales of ~ 250/i -1 Mpc has recently been found by Blanton 
et al.© from the SDSS commissioning data. Similar walls on scales of ~ 250/i _1 
Mpc have been found from the Las Campanas and 2dF redshift surveys near the 
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Northern and Southern Galactic Caps. These results may imply that there 

is a local void of radius ~ 200 - 300/i _1 Mpc in which we live. 

In this paper we analyze directly the SNIa data of the above-mentioned two 
groups using our inhomogeneous models with a local void on scales of about 200 
Mpc, and obtain the confidence contours as well as the best-fit values of the model 
parameters. In §2 we give the basic formulation of the distance modulus and the x 
square for statistical fitting. In §3 we derive the best-fit values of model parameters 
and the confidence contours. In addition, we compare the behavior of our models 
with the new data for z = 1. 7,11a) and we find that they appear to be naturally 
consistent with it. In §4, discussion and concluding remarks are given. 



§2. Distance modulus and the % square 



The theoretical distance modulus is defined by 

/ig = 51og(d L /Mpc) +25, (2-1) 

where g?l is the luminosity distance, related to the angular-diameter distance g?a by 

d L = (1 + z) 2 d A , (2-2) 

along the light ray to a source S with redshift z. This distance modulus is compared 
with the observed one, given by 

Ho = tub - M B , (2-3) 

where Mb and m b are the peak absolute magnitude and the corresponding apparent 
magnitude of a standard SNIa in the B band, respectively. 

In this paper we consider spherical inhomogeneous cosmological models which 
consist of an inner homogeneous region, V , and an outer homogeneous region, V 11 , 
with a boundary of radius ~ 200 Mpc. The observer's position O deviates generally 
from the center C, but is assumed to be close to C, compared with the boundary. 
For the off-center observer, the angular-diameter distance c?a depends not only on z, 
but also on the ajigle between the vectors CO and CS, where S is the source. In our 
previous papers, EP'B 1 the behavior of distances for off-center observers was studied. 
Since the angular dependence is small for remote sources, however, we can ignore it 
for simplicity, assuming that we are approximately at the center. Then d\ depends on 
the source redshift z, the inner and outer Hubble constants Hq and Hq 1 , the inner and 
outer density parameters Q\ and J?q , the outer A parameter Ag 1 , and the boundary 
redshift Z\. The inner A parameter X l is related to Aj, 1 by Aj, = X^(H^/Hl) 2 . The 
equations to be solved to derive d& and the junction conditions are Eqs. (5) - (11) 
in Ref. 10). 

The best-fit values for the cosmological parameters are determined using the 
X-square expression 

X 2 =^2 ( z i l-^o> Hq, ^0) ^o> V) z i) ~~ Mo,-, 
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/(^o,i + ff L,i). ( 2 -4) 

where a^o is the measurement error of the distance modulus and a mz is the dispersion 
in the distance modulus corresponding to the dispersion of the galaxy redshift o~ z 
(coming from peculiar velocity and uncertainty). a mz is related to a z as 

Next, the probability distribution function (PDF) must be obtained to derive 
the confidence contours and the most likely values of the model parameters. It is 
expressed as 

p(Hl Hi 1 , fig, fig 1 , Aq 1 , ^xImo) oc exp ( - ^ , (2-6) 

which is consistent with Eq. (9) in Ref. 2). 

In the present inhomogeneous models there are six parameters, and it is too 
complicated to consider all of them simultaneously. For this reason, we treat here 
several cases with specific values of z\, Hq 1 / Hq, and fig. Then, the corresponding 
normalized PDF can be expressed as 



p(Hl, fig 1 , ASVo) 



ex P (-k 2 ) (2 . 7) 



" dHl JZo d\$ exp(-i X 2 ) ^n 

for given zx,Hq/Hq, and fig, where is the set of distance moduli used and fioz 
is the assumed lower limit of fig 1 , being or — oo. Physically, the region in which 
fig 1 < is meaningless, although it is significant statistically. To obtain a probability 
distribution independent of Hq, we can consider the following quantity: 

/oo 
p(Hln$X\to)dHl (2-8) 
-oo 

The contours are determined by p(fig ! , A^IjUo)- 



§3. Supernova data and the fitting 



At present, we have two SNIa data sets to be used for determining cosmological 
parameters, that of HSST iSchmidt et al.&, Riess et al.&) with 50 SNIa and that 
of SCP (Perlmutter et alJ)) with 60 SNIa. 18 SNIa are common to both data. 
In the former, there are two types of fio data, that due to the Multicolor Light 
Curve Shape method (MLCS) and that due to the template-fitting method. In each 
method, we can use the data for Zj,/io,i and cr M o,i (with i = 1 - 50) given in their 
tables. Following Riess et al.cP for o~ z , we adopt a z = 200 km s _1 and add 2500 
km s _1 in quadrature to a z for SNIa whose redshifts were determined from broad 
features in the SN spectrum. 

In Ref. 4), the data for m^, Zi, <T mz ,i and a Zi i (i = 1 - 60) are given, but the 
values of /io,j were not published and are not available. As shown by Wang,© 
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Fig. 1. The 68.3 % and 95.4 % confidence contours in the " ^o 1 plane in the case 

{zi,Ho/Hl, J?o) = (0.080, 0.82, A). The dotted curves represent the data of Riess et al. (1998), 
the dashed curves those of Perlmutter et al. (1999), and the solid curves the combined data. 



however, we have for 18 SNIa observed by both teams 



AT 



MLCS 



MLCS 
^0 



-19.33 ±0.25, 



(3-1) 



where mff is the effective B-band magnitude of SNIa given by Perlmutter et al., and 
^mlcs - g correS ponding data for due to MLCS method in Riess et al. B 

For the corrected B-band peak absolute magnitude derived by Hamuy et al., 
on the other hand, we have 



M MLCS 



-0.047 ±0.270, 



(3-2) 



which is sufficiently small, compared with the counterpart due to the template-fitting 
method. Accordingly, we adopt the MLCS data of Riess et al.'s two data sets, and 
derive the values of /Uq from the data set of Perlmutter et al. using the relation 



= mf - M B with M B = -19.33, (3-3) 
following Wang.© 

From the two data sets we make a combined data set consisting of 50 SNIa (from 
Riess et al.) and 42 SNIa (from Perlmutter et al.), in which we adopt the data of 
Riess et al. for 18 common data (due to Wang's procedure^). 
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Fig. 2. The 68.3 % and 95.4 % confidence contours in the J7J 1 - A?, 1 plane in the case 
(zuHg/Hl, nl) = (0.080, 0.82, B). 

The above three kinds of data were compared with the theoretical models using 
the following conditions concerning z\, Hq 1 /Hq, and J?q. For zi, we mainly used z\ = 
0.080, corresponding to the radius 240/Zi 1 Mpc, where H\ = lOO/i 1 km s _1 Mpc -1 . 
We compare the PDF in this case with the PDF in cases with z\ = 0.067 (czi/Hq = 
200/Zi 1 Mpc) and z x = 0.100 (czi/flj = 300/Zi 1 Mpc) below. For Hff/Hl, we consider 
three cases, with tig/Hi = 0.80,0.82 and 0.87 (or Hq/H^ = 1.25,1.20 and 1.15), 
respectively. 

For we first consider the inner low-density case (A), in which 

Ql = 0.3 for Q l J > 0.6 (3-4) 

and 

Q\ = for < 0.6. 

For comparison, we also consider the equi-density case (B) with 

Since p J oc (#q) 2 I?q (j = I, II), we have pi = p$ in this case. 

In Figures 1 - 5, we plot the 68.3 % and 95.4 % confidence contours in the 
^q 1 - Aq 1 plane. In Figs. 1, 2 and 3, we use three kinds of data and treat the three 



(3-5) 
(3-6) 
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Fig. 3. The 68.3 % and 95.4 % confidence contours in the J7J 1 - A?, 1 plane in the case 
{zi,Hl l /Hl, Ql) = (0.080, 0.82, A) with fi m = 0. 



cases { Zl ,H^/Hl,Ql) = (0.080, 0.82, A), (0.080, 0.82, B), and (0.080, 0.82, A), with 
J?02 = 0, respectively. Only in the case i?oi = 0, the region of < is excluded. 
It is found comparing these figures that the difference between the two cases A and 
B is small, and the difference between the two cases with l?oz = — oo and is also 
small. 

In Figs. 4 and 5, we display the difference between the contours in the cases 
R = H^/Hl = 0.80, 0.82 and 0.87 and z x = 0.067, 0.080 and 0.100, respectively, 
using the combined data. It is interesting that as R decreases, the contours move 
in the direction of increasing f2ff and decreasing Aq 1 , and for R = 0.80, the flat case 
with vanishing cosmological constant (the Einstein-de Sitter model) is at the center 
of the contours. 

In Table I, the best-fit values of the cosmological parameters in the homogeneous 
models are listed for later comparison. In Tables II, III and IV, we display the best- 
fit values of h 1 , l?^ 1 and Aq 1 with la error bars, and the corresponding values of xt i n 
the cases H^/Hq = 0.82, 0.87 and 0.80, respectively, assuming z\ = 0.080 and Q\ 
(A). The error bars on h l include only statistical errors, not the contributions from 
the errors of the absolute magnitudes. The values of xt fo r the data of Perlmutter et 
al. is found to be much larger than that for the data of Riess et al. This results from 
the situation that for Perlmutter et al. we made artificially the data for /j,q, using the 
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Fig. 4. The 68.3 % and 95.4 % confidence contours in the f?" - Aq 1 plane for the combined data. 
The dotted curves, solid curves and dashed curves represent the cases Hi 1 /Hi — 0.80, 0.82, and 
0.87, respectively. 



Table I. Model parameters determined for the homogeneous models. Here, xl is the value of \ 2 
per degree of freedom. Ho = lOO/i km s _1 Mpc -1 . The errors on h expressed here are only 
statistical. 



data 


Riess et al. 


Perlmutter et al. 


combined 


h 

n 

Ao 

xl 


0.65 ±0.01 
0.2 ±0.6 
0.7 ±0.8 
1.11 


0.666 ± 0.02 
1.0 ±0.4 
1.7 ±0.6 
1.57 


0.65 ±0.01 
0.7 ±0.4 
1.2 ±0.5 
1.45 



Table II. Model parameters determined in the case (21 , Hi 1 / Hi, Q\) = (0.080, 0.80, A), xl is tne 
value of x 2 P er degree of freedom. Hi — WOh 1 km s _1 Mpc -1 . The errors on h 1 are only 
statistical. 



data 


Riess et al. 


Perlmutter et al. 


combined 


ft 1 


0.64 ±0.01 


0.64 ± 0.02 


0.64 ±0.01 


ol 1 


0.1 ±0.7 


0.9 ±0.6 


0.5 ±0.5 


Aj, 1 


-0.9 ± 1.0 


-0.1 ±0.8 


-0.5 ±0.7 


xl 


1.05 


1.62 


1.42 


)fl at 


1.0 ±0.2 


1.0 ±0.2 


1.0 ±0.1 
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Fig. 5. The 68.3 % and 95.4 % confidence contours in the f?" - Aq 1 plane for the combined data. 
The dotted curves, solid curves and dashed curves represent the cases z\ = 0.067, 0.080, and 
0.100, respectively. 

Table III. Model parameters determined in the case (zi,Hq/HI, Q l ) = (0.080, 0.82, A), xl is the 
value of x 2 P er degree of freedom. H\ = lOO/i 1 km s _1 Mpc -1 . The errors on h 1 are only 
statistical. 



data 


Riess et al. 


Perlmutter et al. 


combined 


xl 


0.64 ±0.01 
0.1 ±0.7 
-0.7 ±1.1 
1.05 


0.64 ± 0.02 
0.9 ±0.7 
0.1 ± 1.0 
1.61 


0.64 ±0.01 
0.6 ±0.5 
-0.2 ±0.7 
1.42 


(A) )flat 


0.9 ±0.2 


0.9 ±0.2 


0.9 ±0.1 



average absolute magnitude. In the data of Riess et al., xt in our inhomogeneous 
models is smaller than that in the homogeneous models (Table 1), and so these 
data seem to be fitted better by our inhomogeneous models or the inhomogeneity of 
Hubble constant, though more parameters are included in our case. To the bottom 
of Tables II - V, we add (^J^flat) the density parameter in the case that the outer 
region is spatially flat. 

In Table V we list the best-fit values in the case of fi\ (B), assuming H^/Hq = 
0.82, z\ = 0.080. From a comparison of Tables II and V, we see that the best-fit 
values are not sensitive to In Table VI, moreover, we list the values for the cases 
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with z\ = 0.067 and 0.100, assuming H^/H\ = 0.82 and Q l (A). It is found from 
this table that, because of the smallest xli the fitting for Z\ = 0.080 is better than 
that in the cases with z\ = 0.067 and 0.100. This result may be connected with 
the fact that the boundary radius of the observed local void is likely to be between 
200/Zi 1 Mpc and 300/Zi 1 Mpc. 

Finally,_we compare the behavior of our models with the new supernova data 
for z = 1.7. i3 In Fig. 6 we show the A(m — M) - log z diagram, where A(m — M) is 
the difference between m — M in each model and that in the empty homogeneous 
model (l?o = Ao = 0). Here we have adopted two models with a local void in 
which (flj 1 , A?) = (1.0,0.0) and (0.6,0.0). In both models we use Q l = 0.3, \\ = 
and z\ = 0.067 as the other parameters, to which the behavior of the models is 
not sensitive. The upper, middle and lower curves in each model correspond to 
R (= Hi 1 /Hi) = 0.80, 0.82 and 0.87, respectively. For comparison, we also display 
the behavior of three homogeneous models with (J?o>Ao) = (0.35, 0.65), (0.35, 0.0) 
and (1.0,0.0), and depict in the diagram the binned data (z < 1) and the new data 
(SN 1997ff) for z = 1.7, which are given in Fig. 11 of Ref. 18). It is found from 
Fig. 6 that our models can be consistent with the new data, in contrast to the 
flat homogeneous model with (J?o 5 Ao) = (0.35,0.65). In particular, the model with 
(^q^Aq 1 ) = (1.0,0.0), which is the Einstein-de Sitter model in the external region, 
seems to naturally accord with all data with R = 0.80 and 0.82. In the model with 
(i^o 1 ) Ao 1 ) = (0-6, 0.0), the case R = 0.87 is best among the three cases. If more data 
for z > 1 are provided, we shall be able to distinguish more clearly the behavior of 
the models with a local void and that of the homogeneous models with dominant 
cosmological constant or dark energy. 

§4. Concluding remarks 

In this paper we derived confidence contours and best-fit parameters for inho- 
mogeneous models with a local void using two sets of SNIa data and the combined 
data, and we found that (1) they are very sensitive to the ratio R of the outer Hub- 
ble constant to the inner Hubble constant Hq, (2) the best-fit outer density 
parameter i?^ 1 increases and the cosmological constant parameter Aq 1 decreases as 
functions of R, and (3) (J7q , Aq 1 ) can be (1, 0) for R « 0.8. It is thus found that the 
existence of our local void may solve the puzzling cosmological- constant problem. 

However, we ignored the directional dependence of the [m, z] relation, which may 
be an important factor, especially for nearby SNIa, since magnitudes of SNIa are 
measured by off-center observers. If the data of angular positions of observed SNIa 
are published, it should be possible to take into account the directional dependence, 
and the fitting may be improved. 

The flux averaging proposed by Wangli^P was not carried out here, but it may 
be important when many high-z data with z > 1.0 appear, because they would be 
much affected by the lensing effect. 

From a comparison with the new data for z = 1.7, we found that our models 
are naturally cosnsistent with it, in contrast to the homogeneous models with accel- 
eration due to the cosmological constant or dark energy. For z > 1, the behavior 
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log z 

Fig. 6. The A(m — M) - logz diagram for models with a local void (V) and homogeneous models 
(H). Here A(m — M) = (m — M) in each model — (m — M) in the empty model {Qq = Ao = 0). 
In the two V models, we have upper, middle and lower curves corresponding to R (= H^/Hq) — 
0.80, 0.82 and 0.87, respectively. For comparison the SNIa binned data (z < 1) and the new 
data for SN 1997ff (z = 1.7) are also depicted in the diagram. 



of curves represents the deceleration of the universe, which depends strongly on the 
equation of state of the constituent matter. In our models with a local void, this is 
mainly pressureless matter, while the homogeneous models with the dominant cos- 
mological constant and hypothetical dark energy have negative pressure comparable 
with the mass energy. More data with z > 1.5 is needed in order to distinguish these 
two types of models more clearly and remove the fluctuations due to the lensing 
effect. 
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Table IV. Model parameters determined for the case (z\,Ho/Hl, J?o) = (0.080, 0.87, A), xl is tne 
value of x 2 P er degree of freedom. Hq — lOO/i 1 km s -1 Mpc -1 . The errors on h 1 are only 
statistical. 



data 


Riess et al. 


Perlmutter et al. 


combined 


ft 1 

a" 

X^ 


0.64 ±0.01 
0.2 ± 1.0 
-0.2 ± 1.1 
1.05 


0.65 ± 0.02 
0.7 ±0.7 
0.7 ±1.1 
1.61 


0.64 ±0.01 
0.6 ±0.6 
0.3 ±0.8 
1.42 


(A) )flat 


0.7 ±0.2 


0.7 ±0.2 


0.7 ±0.1 



Table V. Model parameters determined for the case (zx,Hq/HI, flj) = (0.080, 0.82, B). xl is the 
value of x 2 P er degree of freedom. H\ = 100/t 1 km s -1 Mpc -1 . The errors on h 1 are only 
statistical. 



data 


Riess et al. 


Perlmutter et al. 


combined 


h l 

A? 
A? 

X^ 


0.64 ±0.01 
0.0 ±0.7 
-0.8 ±1.1 
1.05 


0.64 ± 0.02 
0.8 ±0.7 
0.0 ± 1.0 
1.62 


0.64 ±0.01 
0.5 ±0.5 
-0.3 ±0.8 
1.42 


(A) )flat 


0.9 ±0.2 


0.9 ±0.2 


0.9 ±0.1 



Table VI. Model parameters determined for the case {z\,Ho/Hl, Q\) = (0.080, 0.82, A), xl is the 
value of x 2 P er degree of freedom. Hq — 100ft 1 km s -1 Mpc -1 . The errors on h 1 are only 
statistical. 



data 


Riess et al. 


Riess et al. 


combined 


combined 


21 
h l 

Mi 

2 

Xv 


0.067 
0.64 ±0.01 
-0.3 ±0.7 
-1.3 ±1.1 

1.10 


0.100 
0.64 ±0.01 
0.4 ±0.7 
-0.3 ±1.1 
1.13 


0.067 
0.64 ±0.01 
0.2 ±0.5 
-0.7 ±0.8 
1.45 


0.100 
0.64 ±0.01 
0.9 ±0.5 
0.2 ±0.8 
1.46 
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